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Abstract 

We study phase transitions in the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) with the weak scale vacuum expectation values of the singlet scalar, con- 
strained by Higgs spectrum and vacuum stability. We find four different types of phase 
transitions, three of which have two-stage nature. In particular, one of the two-stage tran- 
sitions admits strongly first order electroweak phase transition, even with heavy squarks. 
We introduce a tree-level explicit CP violation in the Higgs sector, which does not affect 
the neutron electric dipole moment. In contrast to the MSSM with the CP violation in 
the squark sector, a strongly first order phase transition is not so weakened by this CP 
violation. 
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1 Introduction 



The origin of matter has been a long standing problem in astrophysics and particle physics. 
Several attempts have been made to explain how the matter-antimatter asymmetry was 
generated starting from the symmetric universe, in the early stage of the universe before 
the nucleosynthesis. Among such attempts, electroweak baryogenesis is a scenario which 
is intimately related to physics at our reachfl]. The scenario of electroweak baryogenesis, 
however, requires extension of the minimal standard model, since the minimal standard 
model with the Higgs boson heavier than the present bound does not give rise to strongly 
first-order phase transition, which is required to realize nonequilibrium state, and CP 
violation in the CKM matrix is too small to generate sufficient baryon number. Hence 
one needs extra sources of CP violation and light bosons, which strengthen the electroweak 
phase transition (EWPT). The minimal supersymmetric standard model (MSSM) is one 
of such extensions of the standard model, which can provide strongly first-order EWPT 
and contains many sources of CP violation. The first-order EWPT is expected in the 
MSSM with a light top squark (mi < m t ), which enhances t> 3 -behavior of the effective 
potential, where v represents the expectation value of the Higgs fields. The behavior is 
also enhanced for smaller tan f3, that is, larger top Yukawa coupling. The present bound 
on the lightest Higgs boson restricts the acceptable parameter space of the model, which 
incorporates the strongly first-order EWPT. Further, the first-order EWPT is weakened, 
when the stop sector CP violation, characterized by lm(fiA t e' td ), is large, where 9 is the 
relative phase of the two Higgs doublets[2]. Thus, we are left with rather restricted range 
of the parameter region available for the baryogenesis in the MSSM. 

Another extension of the standard model is the next-to-minimal supersymmetric stan- 
dard model (NMSSM), which contains an extra singlet chiral superfield. It shares the 
advantages of the MSSM and solves the //-problem inherent in the MSSM. The NMSSM 
is reduced to the MSSM in the limit where the singlet decouples leaving a finite //-term. 
For such peculiar parameters, we expect the same behavior of the EWPT as that in the 
MSSM, which has been extensively studied. Here we shall focus on the parameter space 
far from the MSSM, that is, the case where the vacuum expectation value of the singlet 
scalar has the magnitude of the weak scale. One naturally expects such parameters, when 
all the soft masses of the scalars are generated at the same scale where the supersym- 
metry is broken. For such parameters, two of the authors studied Higgs spectrum with 
and without CP violation. In [2], two distinct allowed parameter regions are chosen after 
imposing conditions for the model to yield the sensible eletroweak-breaking vacuum. The 
lightest Higgs boson for one type of the allowed parameter sets has moderate coupling to 
Z boson, while its mass is heavier than 114GeV. The other type of the parameter sets 
contains the Higgs bosons lighter than the bound, but their couplings to the Z boson are 
too small to be produced in the e~e + colliders. For such parameters, only the Higgs boson 
heavier than the bound has the chance to be observed in future collider experiments. 

Our purpose is to study finite-temperature phase transitions in the NMSSM with and 
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without CP violation in the Higgs sector. We find that the transitions in the case with light 
Higgs bosons are very different from that in the others. This paper is organized as follows. 
In Section 2, we introduce the model and define the parameters in the Higgs potential in 
a manner independent of the phase conventions. The effective potential including at the 
one- loop level is given in Section 3. Based on the symmetry of the effective potential, we 
discuss the possible phases and transitions among them in Section 4. In Section 5, we 
present the numerical results on the phase transitions with and without the explicit CP 
violation in the Higgs sector. Section 6 is devoted to discussions on our results and their 
implication on the baryogenesis. We summarize some formulas to calculate the effective 
potential in Appendix. 



2 The model 

The NMSSM is an extension of the MSSM, which contains one singlet superfield N. The 
superpotential of the model is 

W = €ij (y b H d Q^B - y t H l u Q j T + Vl E\VE - \NW d Hi) - (2.1) 

which contains no dimensional parameter. The /z-parameter in the MSSM is induced 
when the scalar component of the singlet N acquires nonzero vacuum expectation value, 
H = A (N). We shall not specify how the supersymmetry is broken, so that we introduce 
generic type of soft supersymmetry-breaking terms. The dimensional parameters such as 
the scalar masses, A-terms and gaugino masses, are considered as inputs, which are con- 
strained, for example, by vacuum stability conditions. The soft supersymmetry-breaking 
terms are composed of those in the MSSM with \x = and the terms including the singlet 
fields: 



*-soft — t'soft - m N n n + 



\A x nH d H u + -A K n 3 + m' N n 2 + h.c. 
o 



(2.2) 



Among these terms the n 2 -term breaks the global Z 3 symmetry, which causes the domain 
wall problem upon broken spontaneously. The n 2 -term is not generated in the simple 
supergravity model, so we shall not include it in the following and consider that the 
discrete symmetry is explicitly broken by some higher dimensional operator at very early 
era before the EWPT. 

We parameterize the expectation values of the Higgs scalars as 

<^> = ( ^ Q Vd ) , <*„> = e ie ( _? Vu "j, (n) = -L e *v n , (2.3) 

where $d, and n are the scalar components of H d) H u and N, respectively. 

The tree-level Higgs potential is written as a function of these expectation values: 



t / 1 22,1 22,1 2 2 

Vo = ■^rn l v d + -m 2 v u + -m N v n 



3 



#2+ #1/2 2\2 _|_ 1^1 ( 2 2. 2 
> V d ~ V u) "+" ■ 



32 



A 



ft 



-e + -e 



2i ¥ >„2 



(2.4) 



where the terms in the first line consist of the soft-supersymmetry breaking terms, while 
those in the second line come from the D- and F-terms. Here we have four complex 
parameters A, ft, A x and A K , and two phases 6q and tpo, which are the values of 9 and tp 
evaluated at zero-temperature vacuum. Some of these are redundant, so that we introduce 
the following notations in order to present the results in a manner independent of phase 
conventions, 



R x = -^Re (\A x e^ + ^) , I x = -J=Im (\A x e^ + ^) , 

Rk = 72 Re ( kA ^ ILP °) ' h = 7f Im { kA ^ ILP °) > 

K = Re (\K*e l ^~ 2 ^) , J = Im (\k* e ue °-^ Q) ) . 



(2.5) 

(2.6) 
(2.7) 



As shown in [2] , these parameters are constrained by the tadpole conditions, which require 
the first derivatives of the potential to vanish at the prescribed vacuum. In particular, 
CP- violation is characterized only by X, since the conditions lead to 



r 1 T T 3 v od v u 3 vl 

2 2 v 0n 2 v 0n 



(2.8) 



where the subscripts denote the values evaluated at the vacuum and tan/3 = Vq u /v 0c1 . 
When all the parameters are real, these equations are reduced to a trivial one. If some of 
them are complex, we must arrange the parameters to respect this condition. In practice, 
we give the absolute values of the parameters as inputs, then choose their phases so as to 
satisfy this condition. As for the CP-even parameters, R x , R K and TZ are related to the 
soft masses. In favor of these parameters, one can eliminate m\, m\ and m 2 N in terms of 
vo n and tan/5o by use of the tadpole conditions. 

The MSSM limit, for which the singlet scalar decouples by taking v n — > oo with Xv n 
and nv n fixed, the model is reduced to the MSSM, in the sense that not only the spectrum 
but also the behavior of the EWPT are the same as those in the MSSM. In the following, 
we focus on the case of weak scale Vo n , for which we expect new features in the Higgs 
spectrum and the EWPT. Then we obtain some new constraints on the parameters in the 
model. In the MSSM, the vacuum parameterized by Vq and tan /3q is the global minimum 
of the potential, as long as the tadpole conditions are satisfied. On the contrary, in 
the NMSSM, we must require that the prescribed vacuum is the global minimum of the 
potential, and masses squared of the other scalars, such as the charged scalars, evaluated 
there are positive. We also impose the condition that the lightest neutral Higgs mass is 
heavier than 114GeV when its coupling to the Z boson is not small. These requirements 
select some allowed parameter sets from the vast parameter space. The allowed parameter 
sets are classified into two kinds. The one contains a light Higgs boson whose mass is 
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smaller than 114GeV and the coupling to Z boson is very small. The other consists of 
the Higgs bosons, all of which are heavier than the bound. We refer to a scenario with 
the former parameter sets as 'light Higgs scenario'. We found in [2] that the allowed 
parameter region of the light Higgs type disappears for tan/3 > 10. As we shall see later, 
the light Higgs scenario is realized for small \k\ and always contain two light bosons. 

In contrast to the MSSM, there is the tree-level CP violation in the Higgs sector, 
characterized by the parameter X. The CP violation in the Higgs sector of the MSSM 
is induced by the loop corrections of the squarks and others. For example, Im (f^A q ) , 
where A q is the squark A-term, generates the scalar-pseudoscalar mixing elements in the 
neutral Higgs mass matrix. This CP violation also affects the neutron EDM, so that 
the squark mass and/or the relative phase of /i and A q are constrained not to exceed 
the upper bound on the EDM. The effect of this CP violation on the Higgs spectrum 
and gauge couplings have been studied by Carena, et al.[3]. They found that sufficiently 
large CP violation mixes the scalar and pseudoscalar, so that the lightest Higgs boson has 
very small gauge coupling to Z boson to escape from the present mass bound. Although 
one expects strongly first-order EWPT with such a light boson, we found that the CP 
violation weakens the EWPT, which was strongly first order because of the light stop in 
the CP-conserving case [4]. 

On the other hand, the CP violation in the Higgs sector of the NMSSM can escape 
from the bound due to the EDM [2], while it can generate the CP-violating bubble wall 
created at the first-order EWPT. We shall see the effect of the tree-level CP violation on 
the EWPT below. 

3 Effective potential 

In addition to the tree- level potential (2.4), we include the one-loop corrections coming 
from the loops of the third generation of quarks and squarks, the gauge bosons and the 
singlet fermion. The mass-squared matrix of the neutral Higgs bosons and the mass of 
the charged Higgs boson are defined as the second derivatives of the effective potential 
evaluated at the vacuum. The explicit forms of these derivaties are given in [2]. 
The finte-temperature effective potential is given by 



V cS (v;T)=V (v)+AV(v;T), 

where V is the tree-level Higgs potential (2.4). The one-loop correction is 

T 4 . (m 7 \ , , , _ T 4 



AV(v;T) = 3 



F„(mJ„) + ^/ 



m, 



+NcE 2E 

q=t,b I j=l,2 



27r 2 "' V T 



m qj 
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(3.1) 
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(3.2) 



where m denotes the field-dependent masses, given in Appendix A and 

og (l - e -^+^) , (3.4) 



F (m 2 ) = 


1 f 

64tt 2 V 


I B (a) = 


/ dx 


Jo 

PCX) 


h{a) = 


/ dx 
Jo 



) • (3-5) 



Here we adopt the DR scheme with the renormalization scale M, which we determined 
in such a way that the one-loop correction to the potential vanishes at the vacuum. The 
order parameters v = (v 1 ,v 2 ,v 3 ,v 4 ,v 5 ) are related to the expectation values of the Higgs 
fields as follows, 

(vi,v 2 ,v 3 ,v 4 ,v 5 ) = (v d , t> u cosA0, w u sinA#, v n cosAip, v n smA(p), (3.6) 
where A9 = 9 — 9 and Ay? — <f — ip . Then the tree- level potential is expressed as 

V (v) = ^mjvl + ^n\{v\ + v\) + ^m 2 N (v 2 4 + v\) 

- [R\{v2Va - v 3 v 5 ) - h(v 3 v 4 + v 2 v 5 )] vi-^ R K {vl - 3vl)v 4 - I K (3vl - vl)v 5 

+ - vl - vl) 2 + ^ [K 2 + vl + vl){vl + vt) + vM + vl 

I |2 



vl + vlf + y K (v 2 (vl - vl) + 2v 3 v 4 v 5 ^j - 1 (v 3 (vj - vl) - 2v 2 v 4 v 5 



(3.7) 



The tadpole conditions recieve the one-loop corrections, which modify the expression 
for the soft masses and the CP-violating parameter I x . The results are summarized in 
Appendix A. 

When we examine the consistency of the parameters and study the EWPT, we search 
for the minimum in the five- dimensional space of the order parameters v. Noting that 
the tree-level potential V and all the field-dependent masses are invariant under 

(^1,^2,^3,^4,^5) !-> (-vi,-v 2 ,-v 3 ,v 4 ,v 5 ), (3.8) 

the effective potential also has this discrete symmetry. Hence, it is sufficient to search for 
the minimum in the space with v u > 0. This symmetry implies that the first derivatives of 
the effective potential with respect to vi, v 2 and 1*3 at the origin vanish at any temperature. 

The phase transitions are studied by searching for minima of the effective potential at 
each temperatur e. Among the phase tra nsitions, the EWPT is defined as the transition at 
which v(T) = ^v\(T) + v 2 (T) + v 2 (T) vanishes when temperature is raised, irrespective 
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of v n (T) = \Jv"l(T) +v%(T), where v(T) denotes the absolute minimum of the effective 
potential at temperature T. For successful electroweak baryogenesis, the EWPT must be 
strongly first order, in order for the generated baryon number not to be washed out by 
the sphaleron process. In the MSSM, this requires a light stop, whose mass is less than 
the top quark mass. Then, the contribution of the light stop to the effective potential, 
proportional to T^I B {fhijT\ produces an effective Tt> 3 -term with a negative coefficient. 
Althought such terms also come from the gauge boson loops, the stop contribution is 
much larger because of the large Yukawa coupling and the color degrees of freedom[5]. 
Such a light stop is realized when the singlet stop soft mass almost vanishes, while the 
doublet soft mass cannot be taken so small at the same time to avoid too large deviation 
of the p-parameter from unity. Hence, a successful baryogenesis in the MSSM requires a 
specific breaking of supersymmetry. 

In the NMSSM, the scalar trilinear term in the tree-level potential is expected to 
yield strongly first order EWPT [6]. Indeed, the EWPT was studied for a wide range of 
parameter space and was found to have more chances to be first order for smaller Higgs 
mass [7]. The naive argument of [6] is as follows: If the EWPT proceeds along the straight 
line in the order parameter space, connecting the origin and the minimum of the effective 
potential corresponding to the broken phase, one can parameterize the order parameters 
as 



v n = y sin a(T), 

with (almost) constant ct(T) and /3(T) at each temperature T. Then the tree-level poten- 
tial is written as 



For appropriate parameters, the coefficient of the y 3 -term is negative, which makes the 
phase transition along y-direction strongly first order. In the MSSM, the EWPT proceeds 
along an almost constant-/? line [8]. Since there is no specific symmetry among the doublets 
and the singlet, the validity of the parametrization (3.9) with a constant a is not obvious. 
As we shall see below, the phase transitions in the NMSSM are classified into several 
types, only one of which admits the parameterization of the order parameters in (3.9). 

4 Possible phases and transitions 

Before presenting the numerical results, we discuss the possible phases and transitions 
among them. A phase at each temperature is characterized by the location of the minimum 



v, 



'u 



v cos P(T) = y cos a(T) cos (3(T), 
v sin (3(T) = ycosa{T) sin/3(T), 



(3.9) 




(3.10) 
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of the effective potential v(T). Among the components of v, v and v n are related to the 
symmetries of the model. Obviously v is the order parameter of the gauge symmetry, while 
v n is the order parameter of a global U(l) symmetry. This is because in the subspace 
of v n — (t> 4 = v 5 = 0), the effective potential is invariant under the global U(l) 
transformation v 2 + iv 3 i— >• e ia (v 2 + iv 3 ). We denote the phases of different symmetries 
as listed in Table 1. For any sensible parameters, the model is in the phase- EW at 



phase 


order parameters 


symmetries 


EW 


V ± 0, V n ^ 


fully broken 


I, I' 


V = 0, V n ^ 


local SU(2) L x U(1) Y 


II 


V ^ 0, V n = 


global U(l) 


SYM 


v = v n = 


SU(2) L x U(1) Y , global 17(1) 



Table 1: Classification of the phases in the NMSSM. 

much lower temperatures than the electroweak scale. At sufficiently high temperatures, 
where all the symmetries are restored, it is in the phase-SYM. In the limit where the 
singlet fields decouple by v n — » oo, the phase- EW at low temperatures transits to the 
phase-I with almost constant v n (T), resulting in the MSSM-like EWPT. In the case of 
Von = O(100)GeV, all the phases can appear at temperature of the electroweak scale so 
that various patterns of phase transitions will be found. 

The existence of the phases-I and II is a novel feature of the NMSSM. In particular, 
in the subspace of v — 0, the potential can develop a nontrivial minimum for some 
parameters. To see this fact, we consider the tree-level potential for v\ — v 2 — v 3 = 0, 

1 1 - |/d 2 

V (v n ) = V (0,0,0,t;„cos Ay^^sin Ay?) = -m 2 N v 2 n - -R K vl + — v n) (4.1) 

where R K = \kA k \ cos[Arg(/tA K ) + 3Atp]. When R 2 K > 4\K,\ 2 m 2 N , V (v n ) has a local 
minimum at v n = a+ = (R K + \J R 2 K — 4 |k| 2 m 2 N ) / (2 |/t| 2 ), where the potential is 

in v <*%( 2 Rl + R^Rl-A\K\ 2 m 2 N \ 

V (a + ) = \m N 12\^ 2 I ' ^ ' 

For a very small \k\ with moderate \A K \ and m 2 N , a + is so large that Vo(®+) becomes 
smaller than the potential at the electroweak vacuum, when the right-hand side of (4.2) is 
negative. That is, for such a parameter set, the phase at zero temperature is the phase-I. 
Although we must exclude such a parameter set, some parameters can admit the phase-I 
of this kind as an intermediate phase at finite temperature. We shall refer to the phase 
realized for a small \k\ as phase-I', in order to distinguish it from the phase-I which occurs 
as the intermediate phase in the MSSM-like phase transition. Both the phases will appear 
at finite temperature with the same symmetry property. However, the Higgs spectrum 
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for the parameters which admits the phase-I' is different from that in the MSSM-like 
parameters, in that the former contains a light scalar. In the limit k — 0, the model is 
invariant under a phase transformation of iV and H u , which is spontaneously broken by 
nonzero v 0n . This symmetry is explicitly broken by small \k\, which results in a light 
Higgs boson. It should be noted that this argument is not so altered by the radiative 
corrections. If we set v± — v 2 — v 3 — in the effective potential, the correction dependent 
on v n arises only from the loop of the singlet fermion, whose contribution is negligible 
for | | <C 1. This also implies that the value of the effective potential in the phase-I' is 
almost independent on temperature. 

As temperature is lowered, the phase-SYM at high temperature transits to the phase- 
EW via various phases. We enumerate the four types of transitions which we encounter 
in the numerical analysis: 

A: SYM -> I =>■ EW B: SYM -> I' =>- EW 

C: SYM II -> EW D: SYM =>• EW 

Here the double arrow indicates the transition at which the electroweak gauge symmetry 
is broken. The second transition of type A is the MSSM-like EWPT, which proceeds 
with almost constant v n (T). The transition of type D was found to be strongly first order 
transition in the case of light Higgs boson[6]. We shall see below that such a transition 
for the parameters consistent with the lower bound on the Higgs mass occurs in the 
presence of a light stop. The other two-stage transitions are novel in the NMSSM. In 
principle, there could be a three-stage transition, but it will be realized for very restricted 
parameters. 

5 Numerical results 

The behavior of the phase transitions is studied by numerically searching for the minimum 
of the effective potential (3.1). We do not use the high-temperature expansion to evaluate 
the integrals (3.5), since there are particles of weak scale mass. The method we employed 
is the same as that in [4]. We perform the minimum search of the effective potential 
for various parameter sets. Among the parameters, the common ones are those known 
from experiments such as Vq = 246GeV, mw = 80.3GeV, mz = 91.2GeV, m t = 175GeV, 
mj = 4.2GeV, and we adopt the A-parameters in the squark sector A t = = 20GeV, 
which are taken not so large that the squark masses- squared do not become negative. As 
for the squark soft masses, we take the following three cases: in the heavy squark case the 
doublet soft mass is taken to be rriq = lOOOGeV and the singlet ones are m~ tR = m~ bR = 
800GeV. In the light squark case-I, rriq = lOOOGeV and m~ tR = m~ bR = lOGeV, and in the 
light squark case-II, rriq = 500GeV and the same singlet soft mass as the case-I. Now the 
remaining parameters are those characterizing the vacuum tan/3 , , ip Q and vo n , and the 
complex ones, A, k, A\, A k . The soft masses in the Higgs potential can be expressed in 
terms of the other parameters by use of the tadpole conditions, as shown in Appendix A. 
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As discussed in Section 3, the independent parameters in the Higgs sector are |A|, \k\, R\, 
R K , 1Z and X. We use the charged Higgs mass m H ±, instead of R\, through the relation 



2 2 1 I a 1 2 2 

m H ± = m w - - \\\ v { 



sin p cos Pq 



1 \ /^AVfoO)" 



2—7^" \ 0#0fc 



(5.1) 



where (• • •) represents the value evaluated at the vacuum. The expression of the second 
derivative with respect to the charged Higgs field is rather lengthy and almost the same as 
that in the MSSM, which is given in [4], except for /x-parameter replaced with \vo n e lipo / V2. 

We thoroughly examined the parameter space of the model for t> 0n — 100 — lOOOGeV in 
the absence of CP violation[2]. The allowed parameter sets are defined to satisfy the two 
conditions; i) the Higgs bosons whose coupling to the Z boson is larger than a tenth of that 
in the standard model Higgs boson be heavier than 114GeV (the spectrum condition), 
and ii) the prescribed vacuum be the absolute minimum of the effective potential (the 
vacuum condition). The second condition also requires that all the masses-squared of the 
scalars in the model be positive, including the squarks and sleptons. The allowed regions 
in the parameter space are roughly classified into two groups; one conains a light Higgs 
boson with small coupling to the Z boson, and the other is similar to the MSSM, with 
the lightest Higgs boson heavier than the present bound. The former exists only when 
Vo n is O(100)GeV, and within the region, \k\ < 0.1. The latter region always exists for a 
large \k\ regime when vo n is small, and extends to small \k\ regime as vo n grows to reach 
the MSSM limit. In the following, we show the results on the phase transitions for several 
points within the allowed regions for v 0n = 200GeV, tan f3 = 5 and A K = — lOOGeV, to 
illustrate the four types of the transitions discussed in the previous section in the absence 
of CP violation. For m H ± = 600GeV in the heavy squark case, the allowed region in 
the (A, k) -plane is shown as a white region in Fig. 1, where the grey region is excluded 
by the spectrum condition and the black region is excluded by the vacuum condition[2]. 
Later, we introduce the explicit CP violation in the Higgs sector which does not induce 
the neutron EDM and examine its effects on the phase transitions. 



5.1 CP-conserving case 

We take all the parameters to be real and set 8q = y?o = 0. For t> 0ra — 200GeV, tan/^o = 5 
and A K = — lOOGeV, we can find all the four types of phase transitions by choosing an 
appropriate set of m#±, A, k and the squark soft masses. The transition of type D is 
observed only with the light squark cases in the light Higgs scenario. The transition of 
type A is found in rather broad region with the heavy Higgs boson, but a light stop 
is necessary for the EWPT to be strongly first order just as in the case of the MSSM. 
The transitions of type B and C appear in the light Higgs scenario at \k\ <C 1. The 
type-C transition also requires a light stop to be strongly first order. Here we take 
(A, k) = (0.9, -0.9) and (0.82, -0.05) with m H ± = 600GeV in the light squark-I case. 
Each corresponds to the transition of type A and C, respectively. An example of the 
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Figure 1: The allowed region in the (A, «)-plane for t> 0n = 200GeV, tan/3 = 5, A K = 
— lOOGeV and m H ± = 600GeV in the heavy squark case. 

type-B transition, we choose (A, k) = (0.85,-0.1) with m H ± = 600GeV in the heavy 
squark case. For an example of the type-D transition, we adopt (A, k) = (0.96, —0.02) 
and rriH± = 700GeV in the case of the light squark-II. We shall refer to each parameter set 
as A, C, B, and D, respectively. The masses and g 2 HZ z °f the mass eigenstates are listed 
inTable 2 for each parameter set. Here guzz is the couplings to the Z boson normalized 
by that of the standard model. For these parameter sets, the only one CP-even scalar has 







Hi 


H 2 


H 3 


#4 


H 5 


A 


m Hi (GeV) 

9H f ZZ 


119.53 


203.59 


265.74 


617.24 


637.47 


0.9992 


5.926 x 10" 4 








1.884 x 10" 4 


B 


m Hi (GeV) 

9H f ZZ 


38.89 


75.31 


131.11 


625.61 


627.95 


6.213 x 10" 8 





0.9999 


6.816 x 10" 5 





C 


m Hi (GeV) 

9H f ZZ 


42.24 


63.49 


117.25 


625.09 


627.44 


0.00188 





0.9980 


9.541 x 10" 5 





D 


m Hi (GeV) 
9uiZZ 


41.88 


58.62.08 


115.15 


730.51 


734.58 





1.015 x 10" 4 


0.9997 


1.632 x 10" 4 






Table 2: The mass and g 2 HZZ of the Higgs mass eigenstates for the four parameter sets. 

the coupling almost equal to unity, while the others have almost zero. This is a general 
feature when the charged scalar is rather heavy, since the tree-level mass matrix of the 
CP-even scalar has an eigenvector near (cos/3o, sin/3o, 0) so that the corresponding mass 
engenstate has g 2 HZZ — 1[2]. The mass eigenvalues are also be understood by recalling 
the properties of the tree-level mass matrices as follows. The mass eigenvalues of the 
pseudoscalar always staisfy m 2 Al < rh 2 < m 2 A ^ where rh 2 = m 2 H± — m\\ + |A| 2 Vq/2. When 
m 2 is large, the mass eigenvalues of the CP-even scalars satisfy m 2 Sl < m| 2 < m 2 < rn 2 Si . 
As discussed in the previous section, there is a light pseudoscalar A\ for small \k\ and 
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m 



7^ 0. Then the relations Tr.M| ~ m 2 (l + v% sin 2 (3 cos 2 l3o/vl n ) 



+ m| and TrTWp 



■Ugsin /9 cos 2 /So/fon) imply that 



m, and that m 5l must be smaller 



than m z , once we impose the condition m S2 > 114GeV. Thus the spectrum of the neutral 
Higgs bosons listed in Table 2 is generic for \k\ ~ 1 (set A) and \k\ <C 1 (B, C, D), when 
the charged scalar is heavy. This should be contrasted to the situation in the MSSM, in 
which the pseudoscalar and the heavier scalar decouple leaving the lighter scalar at the 
weak scale, in the limit of large m#±. 

The phase transitions are studied by searching for the global minimum of the effective 
potential, which is a function of (i>i, v 2, V4) = (vd, v u , v n ) in the CP-conserving case. Which 
type of transitions is realized is anticipated from the structure of the zero-temperature 
potential. To show the potential schematically, we define the reduced effective potential 
as a function of (v, v n ), which is a section of the three-dimensional order parameter space 
with a constant tan/5(T) at the minimum, 



V eS (v, v n ; T) = V eS (v cos (3(T),v sin/3(T), 0, v n , 0; T) - V eS (0, 0, 0, 0, 0; T). 



(5.2) 



The contour plots of the reduced potential at zero temperature for the parameter sets A- 
D are shown in Figs. 2 and 3. The left-hand plot in Fig. 2 suggests that the phase 



A,=0.96, K=— 0.02 



V„ '50 




Figure 2: The contour plot of the reduced effective potential at T = for the parameter 
set A (left-hand side) and D (right-hand side). 

transition proceeds along almost constant v n , while the right-hand plot seems to justify 
the parametrization of (3.9) with a constant a(T). The left-hand plot in Fig. 3 shows there 
is a local minimum along v = 0. As discussed in the previous section, the height of that 
minimum is almost indepenedent of T, while that at the vacuum grows as temperarture 
is raised, so that the transition will be of the type-B. To show the global structures, we 
adopt a larger scale for the vertical axes in Fig. 3 than those in Fig. 2. If we used the 
same scale for them, the contour plot for the set C would looks like that for the set D. 

To see how the transitions proceed, we show the temperature dependence of the local 
minima of the effective potential corresponding to the phases, subtracted by the value at 
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the origin. In the numerical search for the minima, we found a few distinct convergent 
points at the same temperature, which are displayed in the plots as different curves. 
For the set A, the result is plotted in Fig. 4, which shows that the phase-EW changes 
to the phase-I at Tq = 120.47GeV. There the order parameters change from (v, v n ) = 
(106.92GeV, 194.23GeV) to (0, 192.75GeV). This implies the EWPT is first order with 



V eff (GeV 4 ) 
-2.0x1 f 



X=0.9, k=— 0.9 



V eff (OeV 4 ) 
-2.15X10 8 - 



X=0.9, k=— 0.9 




Figure 4: The effective potentials at the local minima corresponding to the phase-EW 
(solid curve) and the phase-I (dashed curve) for the parameter set A. The right-hand plot 
is the close view of the left-hand one near the transition temperature. 

vc/Tc = 0.89. Although we cannot determine the sphaleron decoupling condition in the 
NMSSM, since we have not found the sphaleron solution in the model, that will be similar 
to that in the MSSM for the set A, whose spectrum and phase transition are MSSM-like. 
Then this parameter will not provide sufficiently strong first-order transition, but the 
transition will become stronger for a lighter stop, as in the case of the MSSM. 

An example for the transition of type B is shown in Fig. 5. At Tq = 110.26GeV, the 
low-temperature phase-EW transits to the phase-I', which gives its place to the phase- 
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SYM at T > 500GeV. There the order parameters changes from (v, v n ) = (208.13GeV, 248.85GeV) 
to (0, 599.93GeV). This is a very strongly first-order transition with respect to the elec- 
troweak order parameter. As discussed in the previous section, the slopes of the curves 
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Figure 5: The effective potentials at the local minima corresponding to the phase-EW 
(solid curve), the phase- 1' (dashed curve), the phase-II (dotted curve) and the phase-SYM 
(thick dotted-dashed curve) for the parameter set B. The right-hand plot is the close view 
of the left-hand one near the transition temperature. 

corresponding to the phase-I and I' are smaller than the others. In particular, the potential 
at the phase-I' is amost independent of temperature because of small \k\. 

Fig. 6 illustrates the transition of type C, in which the low-temperature phase-EW 
changes to the phase-SYM through the intermediate phase-II. At Tjv = 98.76GeV, the 
phase-EW changes to the phase-II, which is converted to the phase-SYM at T c = 107.44GeV. 
The order parameters changes from (v,v n ) = (194.27GeV, 173.75GeV) to (165.97GeV, 0) 
at T/v, and from (109.54, 0) to (0, 0) at Tc- The latter is a strongly first-order EWPT. 
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Figure 6: The effective potentials at the local minima corresponding to the phase-EW 
(solid curve), the phase-II (dotted curve) and the phase-SYM (thick dotted-dashed curve) 
for the parameter set C. The right-hand plot is the close view of the left-hand one near 
the transition temperature. 

For the parameter set D, the phase transition proceeds as shown in Fig. 7. Only the 
EWPT occurs at Tc = 103.14GeV where the order parameter changes from (v,v n ) = 
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(182.49GeV, 192.26GeV) to (0, 0). This is a strongly first-order transition first studied in 
[6]. Although there is also the local minimum corresponding to the phase-II in this case, 
it does not take part in the phase transition. 




Figure 7: The effective potentials at the local minima corresponding to the phase-EW 
(solid curve), the phase-II (dotted curve) and the phase-SYM (thick dotted-dashed curve) 
for the parameter set D. The right-hand plot is the close view of the left-hand one near 
the transition temperature. 

The three parameter sets B, C and D are within the allowed region with light Higgs 
bosons. What discriminates the types of the transitions are the appearance of the phase- 
I' and the relative magnitudes of the effective potential at the minima. These factors 
depend on the parameters in the model. All types of the transitions except for the type-B 
requires a light stop to admit the strongly first-order EWPT. This is because the term 
which behaves as Tv 3 with a negative coefficient is effective to make the transition stronger 
for the three types, while the transition of type-B requires a local minimum along v — 
and proceeds by leaping between the minimum and that corresponding to the phase-EW 
so that the u 3 -term is not necessary. 



5.2 CP-violating case 

The phases of the Higgs fields around the phase boundary are important ingredients which 
determine how much baryon number is generated at the first-order EWPT[1]. We found 
that the explicit CP violation in the squark sector weakens the v 3 -behavior of the effective 
potential[4], and expect that the same applies to the transitions which require a light stop. 
Here we study the effect of CP violation in the tree-level Higgs sector characterized by X. 
There are infinite sets of CP violating parameters which yield the same value of X, so we 
constrain them in such a way that the phase relevant to the neutron EDM, ArgA + #o + y?o ; 
vanish 1 [2]. In practice, we set 8q = (po = and take ArgK as an independent parameter, 
from which ArgA K is determined by the conditions (A. 12) and (A. 13). 

1 The EDM is an odd function of this phase plus the phase of the A-teim or the gaugino mass. 
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Here we concentrate on the effect of the CP violation on the transition for the pa- 
rameter sets near the set B, for which the EWPT is strongly first order even with heavy 
stops. For the parameter set B, we introduce the relative phase of k to k = —0.1, 5 K , 
and repeat the numerical search for the minima at each temperature. We find that the 
effective potential at the minimum corresponding to the phase-I' decreases with 5 K and 
becomes smaller than that in the phase-EW even at zero temperature for 5 K > 0.2tt. To 
see this behavior, we adopt another parameter set for which A = 0.83 and k = —0.07 
while the others are the same as the set B, and study temperature dependence of the 
effective potentail at the local minima. The masses and the couplings to the Z boson of 
the Higgs mass eigenstates are listed in Table 3. Since \k\ is not so large, the effect of S K 
on the mixing of the CP-eigenstates is also small. The behaviors of the effective potential 
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m Hi (GeV) 
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43.21 
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66.95 
0.00531 


131.38 
0.9946 


625.40 
6.132 x 10" 5 
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Table 3: The mass and g 2 HZZ of the Higgs mass eigenstates for (A, k) = (0.83, —0.07) and 



5 K = 0, O.lvr and 0.2vr. 



are shown in Fig. 8. For 5 K > 0.37T, the zero-temperature vacuum is in the phase-I'. The 
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Figure 8: The effective potentials at the local minima corresponding to the phase-EW 
(solid curve), the phase-I' (dashed curve), the phase-II (dotted curve) and the phase-SYM 
(thick dotted-dashed curve) for (A, n) = (0.83, —0.07) and 5 K = 0, 0.17T and 0.27T. 

transitions for 5 K = and 0.17T are of type C, while that for 5 K = 0.2tt is of type B. At 
the transition temperatures, the order parameters changes as follows: 



5 K = 0; at T N = 133.22GeV, (v, v n ) = (180.74GeV, 195.49GeV) -> (163.16GeV, 0) 
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at T c = 158.27GeV, (v, v n ) = (26.33GeV, 0) -> (0, 0) 

S K = O.Itt; at = 136.58GeV, (u,u n ) = (173.99GeV, 195.88GeV) -> (154.78GeV, 0) 
at T c = 158.27GeV, (u, u n ) = (26.33GeV, 0) -> (0, 0) 

8 K = 0.2vr; at T c = 120.16GeV, (v, v n ) = (200.62GeV, 208.93GeV) -> (0, 750.93GeV) 

The EWPT of the type B for 5 K = 0.2n is strongly first order as that for the parameter 
set B. This behavior of the potential at the minimum in the phase-I' can be understood as 
follows. As discussed in the previous section, the potential in this phase can be approx- 
imated by the tree-level potential, which is given by (4.2), even at finite temperatures. 
In that expession, both R K and m 2 N are decreasing functions of 8 K < n/2, and V (a + ) 
decreases with S K 2 . 

The phase difference relevant to the baryogenesis is that at the phase transition after 
which the sphaleron process decouples. We have considered such a transition to be the 
EWPT at which v(T) acquires a nonzero value. Even if the EWPT is strong enough to 
suppress the sphaleron process after that, no baryon number can be generated at Tc for 
the transition of type C. As we discussed in the previous section, there is the global £7(1) 
symmetry in the phase-II, so that the phase of V2(T) +ivs(T) is undetermined even in the 
presence of the explicit CP violation. In fact, we observed many degenerate local minima 
with the same \v2(T) +ivs(T)\ but with different phases in the minimal search. This 
implies that the bubble with different phases of the Higgs doublets are randomly created 
at the EWPT. No net baryon number is generated, if such bubbles of a macroscopic 
number are created as suggested by a naive argument of the transition [9]. If v(T) at 
T c is too small to suppress the sphaleron process, just as in the case of 5 K = and 
0.l7r above, one may think that the second transition at Tjv will be another chance for the 
baryogenesis. The same argument, however, applies to this case, so that the relative phase 
of the doublet, as well as that of the singlet, between the two phases are undetermined. 
Thus the transition of type C is not adequate for the baryogenesis. On the other hand, 
the phase of the singlet is determined in the phase-I' in the transitions of type B, since 
the symmetry under phase rotation of the singlet field in the subspace of v i = v 2 = v 3 = 
is explicitly broken by the i? K -term in the tree-level potential. In the example of the 
transition for 5 K = 0.2n above, the <f(T) changes from 0.0492 (phase-EW) to 0.2154 
(phase-I'). Further the phase y?(T) in the phase-I' is independent of temperature, since 
all the field-dependent masses are independent of ip = tan _1 (f 5 /f4) in this phase, so that 
no (^-dependent terms are generated at finite temperatures. 

2 Note that m 2 N is given by (A. 10), while we fixed \k\, |A| and the charged Higgs mass, from which R\ 
is determined. 
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6 Discussions 



We studied various phase transitions in the NMSSM, characterized by the order param- 
eters v(T) and v n (T), for the parameters consistent with the lower bound on the Higgs 
boson, when the expectation value of the singlet scalar takes the value of the weak scale. 
The phase transitions are classified into four types, one of which is the MSSM-like one, 
and the others are peculiar to the NMSSM. These novel three types of transitions are 
realized within the allowed region with a light Higgs boson. One of them is the transition 
studied in [6], but the other two types are equally possible. We found that the transition 
of type B can be strongly first order without need of a light stop. For the parameter sets 
with this type of phase transitions, the Higgs bosons which is expected to be observed 
in future experiments is heavier than the present mass bound and has almost the same 
coupling to the Z boson as the Higgs boson in the standard model. Although the two 
light Higgs boson have very small couplings to the Z boson, their Yukawa couplings to 
the b quark are the same order as the standard model Higgs boson, so that they may be 
observed in hadron collider experiments. In order to obtain the sphaleron decoupling con- 
dition, one must know the sphaleron solution and its energy for the boundary conditions 
corresponding to the relevant phase transition. A work in this direction is in progress and 
the result will be published elsewhere [10]. The results suggest that the sphelaron does 
exist and its energy is a bit smaller than the original sphaleron in the standard model [11], 
because of the negative cubic terms in the potential. 

We also invetigated the effect of the tree-level CP violation on the phase transitions 
of type B and C. Because of the global symmetry in the intermediate phase-II, the model 
which exhibits the transition of type C is not suited for the baryogenesis. On the other 
hand, the difference of <p(T) above and below the transition temperature can be large with- 
out affecting the neutron EDM. This phase becomes dependent on the spatial coordinate 
near the bubble wall created at the first-order transition. Noting that the /z-parameter 
in the MSSM is effectively induced as /j, — \v n e lip /y/2, this introduces a new source of 
space-dependent phase into the mass matrix of the charginos, neutralinos, squarks and 
sleptons, which will produce new contributions to the source of the baryon number in the 
super symmetric standard models[12]. 

In this work, we have been focused on the static properties of the phase transitions. To 
evaluate the generated baryon number, we must know about their dynamics. Although 
the study of the dynamical aspects of the transitions is beyond the scope of this work, 
we expect that as long as all the mass parameters in the Higgs potential are of the order 
of the weak scale, the dynamics is not so far from those in the standard model or the 
MSSM. In fact, we have not encountered an extremely deep minimum or high barriar 
between the minima of the effective potential in the numerical studies, which may delay 
the transition or eternally trap the model within an unphysical state. This situation might 
be altered, if we consider the case of larger v 0n but not so large as in the MSSM limit. It 
will be interesting to investigate the phase diagram of the model for broader range of the 
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parameters. 



Acknowledgements 

The authors gratefully thank A. Kakuto and S. Otsuki for valuable discussions. This 
work was supported by the kakenhi of the MEXT, Japan, No. 13135222. 



A Notations 

A.l field-dependent masses 



Here we summarize the field-dependent masses m used to define the effective potential in 
terms of the vector v defined by (3.6). Those of the standard model particles are written 
by the components in the doublets as 



-2 02 + 01 i 2 i 2 i 2\ 

m? = (uj +v 2 +v 3 ), 



™2 



'■Z 
.2 



m 



w 



(A.l) 

The those of the squarks and the singlet fermion depends on the singlet component (t> 4 , v 5) 
and are given by 



m. 



\Vt / 2 1 2\ 

(V2 + V 3 ), 
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- 2 o I I 2 / 2 1 2\ 

m ^ N = 2|K| (U4+V5), 



(A.5) 



where 



1 / 2 5 o 



(A.6) 



R q = l=Re(XA q e^ + ^), I q = -±=1m (XA q e^ + ^) , (q = t,b). (A.7) 
As for the convention for the squark sector, refer to [4]. 
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A. 2 one-loop tadpole conditions 



The mass parameters in the tree-level potential is written in terms of the others by 
requiring that the first derivatives of the zero-temperature effective potential evaluated 
at the vacuum be vanish. We call these conditions as the tadpole conditions. The three 
conditions derived from the derivetives with respect to the CP-even fields are used to 
express mf, m\ and m 2 N as 
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The first derivatives with respect to the CP-odd fields lead to the relation among the 
CP- violating parameters as 
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(A.12) 
(A.13) 



Here ^ and 6j are the components of the vector i and 6 defined in [2], respectively, and 
Umlml) = m 2 (log^| - l) ± (log^| - l) . (A.14) 
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